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Abstract 

We present the results of a search for high energy extraterrestrial neutrinos with 
the Baikal underwater Cherenkov detector NT200, based on data taken in 1998 - 
2003. Upper limits on the diffuse fluxes of v e + + z/ r , predicted by several models 
of AGN-like neutrino sources, are derived. For an E~ 2 behavior of the neutrino 
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spectrum, our limit is E 2 <& U (E) < 8.1 x 10 7 cm 2 s 1 sr 1 GeV over an neutrino 
energy range 2 x 10 4 -f 5 x 10 7 GeV. The upper limit on the resonant v e diffuse flux 
is $ Pe <3.3xl(T 20 cm-W^Gey- 1 . 

Key words: Neutrino telescopes; Neutrino astronomy; UHE neutrinos; BAIKAL 
PACS: 95.55.Vj; 95.85.Ry; 96.40.Tv 



1 Introduction 

High energy neutrinos are likely pro- 
duced in many violent processes in 
the Universe. Their detection would 
unambiguously reveal the hadronic 
nature of the underlying processes. 
Neutrinos would be generated by 
proton-proton or proton-photon in- 
teractions followed by production 
and decay of charged mesons. 

One detection mode of neutrino tele- 
scopes is the identification of individ- 
ual, point-like sources of high energy 
neutrinos. Galactic candidates for 
these objects include supernova rem- 
nants or microquasars, extragalactic 
objects are e.g. active galactic nu- 
clei (AGN) and Gamma Ray Bursts 
(GRB). Individual sources might be 
too weak to produce an unambiguous 
directional signal, however the inte- 
grated neutrino flux from all sources 
could produce a detectable diffuse 
neutrino signal. Astrophysical neu- 
trinos generated in top-down models 
are, by definition, of diffuse nature. 
A diffuse neutrino flux can be iden- 
tified by a high-energy excess on top 
of the background of known fluxes of 
charged particles recorded by a neu- 
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trino telescope. Such charged parti- 
cles are dominantly muons produced 
in the atmosphere above the detec- 
tor, with a small contribution from 
muons generated in interactions of 
atmospheric neutrinos. 

In this paper we present results of 
a search for neutrinos with energies 
larger than 10 TeV. The analysis is 
based on data taken with the Baikal 
neutrino telescope NT200 between 
April 1998 and February 2003. In- 
stead of focusing to high-energy par- 
ticles crossing the array, the analysis 
is tailored to signatures of isolated 
high-energy cascades in a large vol- 
ume around the detector [1]. This 
search strategy dramatically en- 
hances the sensitivity of NT200 to 
diffuse high energy processes. 

The cascades can stem from leptons 
and hadrons produced in high energy 
charged current processes 

ui(ui) + N ^ /-(/+) + hadrons, (1) 

or from the hadronic vertex of neutral 
current processes 

NC 

vi{vi) + N — > vi(yi) + hadrons, (2) 

where I — e, // or r. The energy re- 
leased by the hadronic cascade in re- 
action (2) is small compared to that 
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of the leptonic cascades in (1). Since 
only electrons and taus develop cas- 
cades (electrons by directly shower- 
ing up, taus via their decay to sec- 
ondary particles which develop a cas- 
cade), the sensitivity of this search is 
dominated by u e and u T detection. 

Cascades can also be produced by 
resonant W production 

z/ e + e~ — > W~ — > anything, (3) 

with the resonant neutrino energy 
E = Ml/2m e = 6.3 x 10 6 GeV and 
cross section 5.02 x lCT 31 cm 2 . 

This paper is organized as follows: 
In section 2 we describe detector and 
site, and illustrate the detector per- 
formance. In section 3 we first line 
out the search strategy. This is fol- 
lowed by a description of how a laser 
light source mimicking high energy 
cascades is used to investigate the 
detector performance with respect 
to energetic cascades. Furthermore, 
we describe the simulation of signal 
and background processes. Section 4 
is devoted to the analysis of the ex- 
perimental data. Limits on high en- 
ergy neutrino fluxes and comparison 
to models are presented in section 
5. Section 6 gives conclusions and 
a short outlook on how to further 
improve the present limits. 



of the telescope, NT200 [2,3,4], takes 
data since April 6th, 1998 and con- 
sists of 192 optical modules (OMs). 
A schematic view of the Baikal Tele- 
scope NT200 is shown in fig. 1. An 
umbrella-like frame carries 8 strings, 
each with 24 pairwise arranged OMs. 
All OMs face downward, with the 
exception of the second and eleventh 
pairs on each string which face up- 
ward. Three underwater electrical 
cables connect the detector with the 
shore station. Each OM contains a 
37-cm diameter QUASAR - photo 
multiplier (PM), which has been 
developed specially for our project 
[5]. The PMs record the Cherenkov 
light produced by charged particles 
in water. The two PMs of a pair are 
switched in coincidence in order to 
suppress background from biolumi- 
nescence and PM noise. A pair of 
OMs defines a channel. The light ar- 
rival time assigned to a channel is the 
response time of the OM with the 
earliest hit. The amplitude assigned 
to a channel is that recorded by one 
pre-selected PM of the two PMs in 
a pair. For those channels for which 
one of the two OMs failed and the 
single PM noise rates of the remain- 
ing OM are not too high (smaller 
than the average noise rate of 100 
kHz), that OM is operated in single 
mode (lPM/channel), with a thresh- 
old slightly higher than the typical 
0.3 photo electrons. 



2 Detector and site 

The Baikal Neutrino Telescope is op- 
erated in Lake Baikal, Siberia, at a 
depth of 1.1 km. The present stage 



A trigger is formed by the require- 
ment of > iV hits (with hit referring 
to a channel) within 500 ns. N is typ- 
ically set to 3 or 4. For these events, 
amplitude and time of all fired chan- 
nels are digitized and sent to shore. 
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To Shore 




Fig. 1. A schematic view of NT200. The expansion left-hand shows 2 pairs of optical 
modules ("svjaska") with the svjaska electronics module, which houses part of the 
readout and control electronics. 
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Fig. 2. Left panel: the normalized amplitude distribution for channel 2 (points - 
experiment, histogram - atmospheric muon simulation). Right panel: the normal- 
ized time difference of channel 42 and channel 43. Histogram - experiment, dashed 
histogram - expectation from atmospheric muons. 



Two nitrogen lasers are used for cali- 
bration of the detector. The first one 
{fiber laser) is mounted just above 
the array. Its light is guided via op- 
tical fibers of equal length to each 
OM pair. The fiber laser provides the 
OMs with simultaneous light signals 
in order to determine the offset for 



each channel. The second laser (wa- 
ter laser) is arranged 20 m below the 
array. Its light propagates through 
the water. This laser serves to mon- 
itor the water quality, in addition 
to dedicated environmental devices 
located along a separate string. In 
the context of this analysis, however, 
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its main purpose is to simulate high 
energy particle cascades outside the 
geometrical volume of NT200. The 
maximum light intensity of the water 
laser (~10 n photons/pulse) roughly 
corresponds to the Cherenkov radia- 
tion emitted by high-energy cascades 
with E cas «1 PeV. A full cycle of 
detector calibration running both 
lasers over a wide range of intensities 
is repeated every third day. 

Figure 2 demonstrates the level to 
which basic features of the detector 
are understood and can be repro- 
duced by Monte-Carlo (MC) simu- 
lations. We have used atmospheric 
muons as a high statistics standard 
calibration signal. The figure shows 
distributions of recorded amplitudes 
and light arrival time differences 
compared to MC-simulation. The 
experimental data are in good agree- 
ment with simulation (see also [2]). 

Lake Baikal deep water is charac- 
terized by an absorption length of 
L abs (480 nm) =20-^24 m, a scatter- 
ing length of L s =30-h70 m and a 
strongly anisotropic scattering func- 
tion f(9) with a mean cosine of the 
scattering angle cos(#) = 0.85 -r- 0.9. 



3 The analysis method 

3.1 Search strategy 

The BAIKAL survey for high en- 
ergy neutrinos searches for bright 
cascades produced at the neutrino 
interaction vertex in a large volume 
around and below the neutrino tele- 
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Fig. 3. Detection principle for neutrino 
induced cascades with NT200. 

scope (see fig. 3). Lack of significant 
light scattering allows to monitor a 
volume exceeding the geometrical 
volume by more than an order of 
magnitude. This results in sensitivi- 
ties for high energy cascade detection 
which are close to those of the much 
larger AMANDA detector [6]. The 
main background source to this anal- 
ysis are atmospheric muons, with a 
flux 10 6 times higher than that of 
atmospheric neutrinos. 

We select events with high multiplic- 
ity of hit channels N^, correspond- 
ing to bright cascades. The main re- 
maining background to isolated cas- 
cades from neutrino interactions are 
then cascades from bremsstrahlung 
along energetic downward muons. To 
separate high-energy neutrino events 
from background events a cut on the 
variable £ mm = min(tj — tj) (with 
i < j) is applied. Here, t i; tj are the 
arrival times at channels i, j on each 
string (the numbering of channels is 
from top to bottom along strings), 
the minimum over all strings is cal- 
culated. Positive and negative values 
of t m j n correspond to upward and 
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downward propagation of light, re- 
spectively. We require 

*min > -10 ns- ( 4 ) 



This cut accepts only time patterns 
corresponding to upward traveling 
light signals. It rejects most events 
from brems-cascades produced by 
downward going muons since the 
majority of muons is close to the ver- 
tical; they would cross the detector 
or pass nearby and generate a down- 
ward time pattern. Only few muons 
with large zenith angles may escape 
this cut and illuminate the array 
by their own Cherenkov radiation 
or that from bright cascades from 
below. 

The energy spectrum of neutri- 
nos from galactic and cosmological 
sources or from the decay of topo- 
logical defects is expected to have 
a significantly flatter shape than 
the spectrum of atmospheric muons. 
This is reflected by different iV^ dis- 
tributions. In fig. 4 we show dis- 
tributions of simulated events which 
survive cut (4) and would be induced 
by electron neutrinos with fluxes of 
shape ~ E^ 1 , with 7=1.5, 2 and 2.5 
(normalized to each other, simula- 
tions is for 80 operating channels). 
Also shown is the distribution 
of background events induced by 
atmospheric muons (normalized ar- 
bitrarily to neutrino events). This 
distribution is much steeper, so that 
an extraterrestrial neutrino signal 
would appear as an excess of events 
with large iV^t above the back- 
ground of atmospheric muons. 
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Fig. 4. The normalized hit multiplicity 
distributions of events from u e fluxes af- 
ter selection criterion (4). Solid, dashed 
and dotted histograms correspond to 
7=1.5, 2, 2.5, respectively. Also shown 
is the iV^u distribution of background 
events from atmospheric muons (thick 
line) . 

3.2 Laser calibration 



As mentioned above, a central pur- 
pose of the water laser is to investi- 
gate the response of NT200 to bright 
light sources below the detector. The 
water laser is mounted on the cen- 
tral string about 20 m below the 
bottom OMs of NT200. An atten- 
uator allows to operate the device 
with five gradually decreasing light 
pulse intensities. In the most power- 
ful mode, about 10 11 photons/pulse 
are emitted isotropically. Figure 5 
(left panel) shows the amplitudes on 
channel 49 (turned downward) and 
channel 50 (turned upward) induced 
by these pulses. The distances be- 
tween the water laser and channels 
49 and 50 (both on central string 
and above water laser) are 88.8 m 
and 81.3 m, respectively. In contrast 
to channel 49, which faces the light 
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Fig. 5. Laser calibration: left panel - the amplitude distributions of channels di- 
rected downward (channel 49) and upward (channel 50). Right panel: time differ- 
ence distribution St for pairs of channels facing both downward (49,51) and upward 
(50)/downward (51), respectively. 



source, channel 50 is directed in the 
opposite direction. Therefore it is 
almost blind to direct light and de- 
tects mostly scattered photons. Due 
to the low scattering coefficient of 
Baikal water, the signals in channel 
50 are about 70 times lower than in 
channel 49. 

Figure 5 (right panel) shows the dif- 
ference St of laser light arrival times 
on channel 49 (channel 50) and chan- 
nel 51. Like channel 49, channel 51 is 
turned downward. It is located 12.5 
m below channel 49. The r.m.s. of 
5t = (t 4 g —t 5 i), i.e. for the two down- 
ward looking channels, is about 1.4 
ns. It is not affected by light scat- 
tering in water and entirely due to 
experimental time resolution. The 
r.m.s. of 5t = (t 50 — t 51 ) is about 8.4 
ns. This large value is caused by the 
dispersion of photon arrival times on 
channel 50 due to light scattering. 
To reach maximum background sup- 
pression in the cascade search by ap- 
plication of strict arrival time cuts, 



we exclude upward looking channels 
from the analysis. 

The distributions of the variables 
N-^ft and t m i n , which are used for the 
selection of neutrino induced events, 
are shown in fig. 6. Distributions on 
the left panel correspond to the five 
laser intensities - 1.2 x 10 11 , 2.2 xlO 10 , 
4.7xl0 9 , 1.0 xlO 9 and 2.5 x 10 s pho- 
tons/pulse - which approximately 
correspond to cascade energies 1200, 
220, 47, 10 and 2.5 TeV, respectively. 
The t m [ n distribution on the right 
panel corresponds to the most pow- 
erful laser intensity mode. In agree- 
ment with above discussion of the 
upward moving light selection crite- 
ria (4) we found for nearly all laser 
events t m i n > —10 ns. 

The high transparency and the low 
scattering of Baikal water allow the 
reconstruction of coordinates and 
intensity of bright light flashes. The 
laser experiment allowed to verify the 
reconstruction procedure of shower 
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Fig. 6. Left panel: the hit channels multiplicity of events induced by laser pulses 
with five different intensities for. Right panel: the i m in distribution for laser events. 



position and intensity. We found a 
vertical laser coordinate precision of 
~1 m and relative intensity precision 
of -10 % [7]. 



3.3 High-energy neutrino simula- 
tion 



The number of expected events N Vi 
during observation time T is: 

N Vi =T J dVtJ dE h V eff (Q,E sh )x 
E / N APH2 o^$ Ui {ti,E u ,X)dE„{5) 

k h 



transferred to the hadron, E s ^ - 
total energy of secondary showers, 
V e g(fl, E^) - the detection volume. 
The index z/j indicates the neutrino 
type 0; = v e , z/ e , Vp, v T , v T ) and 
k =1,2 corresponds to CC- and NC- 
interactions respectively. N A is the 
Avogadro number and Ph 2 o the wa- 
ter density. 

For v e and u^, the flux <& v satisfies the 
following transport equation: 



d$ v {E v 



dX 

dKMK 



-N A <T tot <£> v (E, 
da NC 



dE„ 



(K,E U ) (6) 



L 

X(ti) = J Pearth{l)dl, 


where <& Vi (Q,E u ,X) is the flux of 
high energy neutrinos with energy 
E v in the vicinity of the detector, 
Q - the neutrino direction, X(Q) 
- the thickness of matter encoun- 
tered by the neutrino on its passage 
through the Earth, E^ - the energy 



where a tot = a cc + a NC . 

For u T , the tau neutrino and tau lep- 
ton fluxes satisfy the equations: 
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flux ratio for neutrino and antinu- 
trino of vjv =1 1 . 



dX \i ec (E T ) 

dE Vr ^{E Vr ) — {E Vr ,E T ) 

where Xf ec = ^cr T p is the decay 
length of the r-lepton, dn/dE Ur - v T 
production probability. The trans- 
port equation for the r-lepton is 
valid under the assumption that 
\f ec / \ T « 1, where A r - the interac- 
tion length of r propagation through 
the medium. This assumption is 
valid for energies E T <10 8 GeV. 

A MC-code is used to solve equa- 
tions (6) - (7), with the bound- 
ary conditions for neutrino fluxes 
$ Vi (E, 0) = A Vi f Vi (E), where f Vi (E) 
is a diffuse AGN-like flux or other 
predicted UHE neutrino fluxes, and 
a A v . a normalization coefficient. For 
tau leptons $ T (E, 0) = 0. For neu- 
trino interactions and tau-neutrino 
regeneration we used cross-sections 
from [8,9]. The neutrinos are prop- 
agated through the Earth assum- 
ing the density profile of the Pre- 
liminary Reference Earth Model 
[10]. Although a flavor ratio of 
v e : : u T PS 1:2:0 is predicted 
for generic neutrino fluxes at cos- 
mic sources, equal fractions of all 
three neutrino flavors are expected 
at Earth because of neutrino oscil- 
lations. Throughout this paper we 
assumed a neutrino flavor ratio at 
Earth of v e : : v T =1:1:1 and the 
same shape of energy spectra f(E u ) 
for all neutrino flavors, as well as a 



The detector response to Cherenkov 
radiation of high energy cascades was 
simulated taking into account the 
effects of absorption and scattering 
of light as well as light velocity dis- 
persion in water [11,12,13]. We also 
implemented the longitudinal de- 
velopment of cascades. For electron 
cascades with i? sn >2xl0 7 GeV and 
for hadronic cascades with E s ^ >10 9 
GeV, the increase in cascade length 
due to the LPM effect [14] was ap- 
proximated as E 1 ^ according to Ref. 
[15]. 

3.4 Atmospheric muon simulation 

Downward going atmospheric muons 
are the most important source of 
background. The simulation chain 
of these muons starts with cosmic 
ray air shower generation using the 
CORSIKA program [16] with the 
QGSJET [17] interaction model and 
the primary composition and spec- 
tral slopes for individual elements 
taken from [18]. Atmospheric muons 
are propagated through the water 
using the MUM program [19]. Dur- 
ing passage through the detection 
volume the detector response to 
Cherenkov light from all muon en- 
ergy loss processes is simulated. 

For illustration of the consistency 

1 A violation of this assumption (e.g. 
for neutrino production in p^ interac- 
tions) has a small influence on the result 
due to the similarity of u and v cross- 
sections in our energy range. 
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Fig. 7. Left panel: the £ m j n distributions for experiment (dots) and atmospheric 
muon simulation (bold histogram) for iVu:+ >40. Histograms 1 and 2 correspond 
to atmospheric muon simulation with cos# > 0.8 and cos 8 < 0.8, respectively. 
Right panel: the distributions for experiment (dots) and atmospheric muon 

simulation (histogram) with i m j n >-100 ns. 



of experimental data with simula- 
tion we present in fig. 7 (left panel) 
the t m [ n distributions for simulated 
background events (histograms), as 
well as for experimental events with 
a hit multiplicity > 40 taken 

during 41.8 live days in 1999. Curves 
1 and 2 correspond to atmospheric 
muons with zenith angles cos# > 0.8 
and cos 6 < 0.8, respectively. We find 
good agreement for all t m i n values, 
with a small deficit only for vertical 
downward going muons in the inter- 
val -300, -230 ns (i.e. far-off the cut 
value -10 ns) . This deficit is due to the 
difficulty to precisely determine the 
low OM detection efficiency of down- 
ward looking OMs for straight back- 
ward light illumination by vertical 
atmospheric muons. Figure 7 (right 
panel) shows the distribution of 
events with £ m j n >-100 ns as well as 
background simulation (histogram). 
We conclude that experimental t m j n 
and A^j t distributions are consis- 



tent with distributions expected for 
atmospheric muons. 



4 Data selection and analysis 

Within the 1038 days of the detec- 
tor live time between April 1998 and 
February 2003, 3.45 x 10 8 events 
with A^j t > 4 have been recorded. 
For this analysis we used 22597 
events with hit channel multiplicity 
^hit > ^ w hich obey the condition 
(4). For these events upward looking 
channels as well as channels which 
are operated in the 1-PM/channel 
mode (see section 2) have been ex- 
cluded from the following analysis. 

During the 1038 days NT200 took 
data in various configurations, which 
have evolved due to groups of OMs 
failing. Neglecting few-OM differ- 
ences, the data can be grouped ac- 
cording to three basic configurations. 
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Fig. 8. Left panel: the t m \„ distribution of experimental events (dots) which survive 
cut (4) as well as the expected distribution of simulated background events (his- 
togram). Right panel: the -/VW distribution of experimental events (dots) as well as 
the background prediction (histogram) after cut (4). 



Table 1 

NT200 effective configurations. 



Conf. 


N ov 


T 

(days) 


AW 


^rnax 


1 


71 


316 


12146 


47 


2 


59 


612 


9473 


42 


3 


46 


110 


978 


32 



The average number of working 
channels N p, the data taking time 
T, the number of detected events 
N eY , as well as the largest multiplic- 
ity of hit channels for detected events 
A^hrf X are shown in Table 1. 

Figure 8 shows the t m [ n and iV^t dis- 
tributions for experiment (dots) and 
background simulation (histograms). 
We conclude that after application 
cut (4) data are consistent with sim- 
ulated background for both t m - m and 
Nfoft distributions. No statistically 
significant excess above the back- 
ground from atmospheric muons has 
been observed. 



To maximize the sensitivity to a neu- 
trino signal we introduce a cut in 
the (t m i n , N-^ft) phase space. Figure 
9 (left panel) shows the population 
of the (i m j n , iVjjit) phase space for 
experimental events (triangles) as 
well as for background simulation 
(dots). The distribution for neutrino 
induced events (dots) is shown in 
the right panel of fig. 9 (for spectral 
index 7 = —2). We note that back- 
ground events populate the lower-left 
part of the plot, in contrast to signal 
events. The bounds which fence the 
(*min> -^hit) area w hich we assigned 
to background are listed in Table 2 
for the three NT200 configurations. 

With no experimental events outside 
the area populated by background 



events in 



the (t 



mini 



iVi 



hit; 



phase 

space we can derive upper limits on 
the fluxes of high energy neutrinos 
as predicted by different models of 
neutrino sources. 

The neutrino detection energy range 
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Fig. 9. Left panel: distributions of experimental (triangles) and expected background 
(dots) events in the (£ m i n , iVj 1 ^)-plane. Right panel: same distribution for events 
expected from high energy cascades. Results are shown for configuration 1 (Table 
1). The line shows cut condition (see Table 2). 

Table 2 



Cut conditions in the (imin, JV^J-plane used for the selection of neutrino induced 
events (train in ns). 



Conf. 1 


-10< t min <5 


5 < t min <30 


*min >3 ° 




N Ut > 50 


N hit > lO^ 1 - 8-0 - 02 *™") 


N Ut > 16 


Conf. 2 


-10 < t mhl <10 


10< t min <30 


*min >3 ° 




N hrt > 44 


N hit > 10( L865_0 - 0223tm ™) 


N hrt > 16 


Conf. 3 


-10 < t min <10 


10 < t miQ <30 


*min >30 




N hit * 33 


iV hit > 42-0.9t min 


*hit > 16 



of NT200 which contains, for in- 
stance, 90% of expected events, de- 
pends on the energy shape of the 
neutrino flux. Figure 10 shows the 
fractions of expected events induced 
by diffuse neutrino fluxes following 
an E' 1 shape with spectral indices 
7=1, 1.5, 2, 2.5 and 3.7, for ener- 
gies ranges 10 TeV< E < E v . The 
effective detection range is shifting 
towards higher energies with de- 
creasing 7. Table 3 shows the energy 
ranges for effective detection of neu- 



trino fluxes with different 7, as well 
as the median values of the energy 
distributions. 

With increasing energy, neutrinos are 
stronger absorbed when propagating 
through the Earth, and upward go- 
ing fluxes are suppressed. Our search 
strategy accepts neutrinos from all 
directions. With increasing neutrino 
energy the angular distributions are 
shifted towards smaller zenith angles. 
Zenith angle ranges A6 which con- 
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Fig. 10. The fraction of expected events 
induced by diffuse v e fluxes with spec- 
tral indices 7=1, 1.5, 2, 2.5 and 3.7 
within energy range 10 TeV< E < E v . 

Table 3 

Detection energy range and median en- 
ergy for different spectra. 

7 AE, TeV £? med , TeV 



1.0 


5xl0 3 - 7xl0 5 


2.2xl0 5 


1.5 


2xl0 2 - 5.6xl0 5 


2.0xl0 4 


2.0 


22 - 5.0xl0 4 


5.6xl0 2 


2.5 


14 - 2.0xl0 3 


63 


3.7 


10 - 89 


20 



tain 90% of expected events induced 
by an E~ 2 neutrino flux are 75° - 150° 
and 40° - 90° for energy ranges 10 
TeV - 100 TeV and 10 5 - 10 6 TeV, re- 
spectively. 



The detection volumes V e fj for all 
three neutrino flavors after all cuts 
were calculated SIS db function of neu- 
trino energy and zenith angle 9. The 
energy dependence of the detection 
volumes, averaged over all neutrino 
arrival directions, are shown in fig. 
11 (left panel). The value of V e Q 



rises from ~10 5 m 3 for 10 TeV up 
to (4-6) xlO 6 m 3 for 10 4 TeV and 
significantly exceeds the geometrical 
volume V g « 10 5 m 3 of NT 200. This 
is due to the low light scattering and 
the nearly not distorted light fronts 
from Cherenkov waves originating far 
outside the geometrical volume. In 
the case of u e detection, the volume 
saturates for E Ve >10 4 TeV because 
of the LPM effect. Figure 11 (right 
panel) illustrates the difference be- 
tween V e -g and Vg. Shown here are 
the coordinates of neutrino interac- 
tion vertices for events which survive 
cuts (4) (dots) and all cuts (rectan- 
gles), assuming an E~ 2 spectrum. 



Systematic uncertainties in the op- 
tical properties of water, the abso- 
lute detector sensitivity and the neu- 
trino cross sections at high energies 
influence the number of expected sig- 
nal events. The uncertainty of 10% in 
absorption length results in 20% un- 
certainty in the number of expected 
events. The uncertainty of 14% in the 
sensitivity of OMs results in 10% un- 
certainty in the number of expected 
events. Below 10 16 eV, all standard 
sets of parton distributions yield very 
similar cross sections. Above this en- 
ergy, the cross sections are sensitive 
to assumptions made about the be- 
havior for x — > 0. The uncertainties 
of cross sections are less than 10% 
below 10 18 eV [8,20,21] and result in 
uncertainties of 8% in the number of 
expected events. Treating these er- 
rors as independent and adding them 
quadratically, the signal uncertainty 
becomes 24%. 
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Fig. 11. Left panel: energy dependence of detection volumes. Right panel: coordinates 
of v e interaction vertices for events which fulfill cut (4) (dots) and all final cuts for 
neutrino event selection (rectangles), respectively. The ordinate gives the vertical 
position, the abscissa the radial distance, with (z=0, R=0) corresponding to the 
center of NT200. The small rectangle marks the geometrical volume of NT200. 

Table 4 



Expected number of events -^ moc j e i and model rejection factors for models of astro- 
physical neutrino sources. The assumed upper limit on the number of signal events 
with all uncertainties incorporated is n go % = 2.5 



BAIKAL 


AMANDA [6,24] 


Model 




«/„ 


v T 




re 90%/^model 


"90%/^model 


1(T 6 x E~ 2 


1.33 


0.63 


1.12 


3.08 


0.81 


0.86 


SS Quasar [25] 


4.16 


2.13 


3.71 


10.00 


0.25 


0.21 


SP u [26] 


17.93 


7.82 


14.43 


40.18 


0.062 


0.054 


SP 1 [26] 


3.14 


1.24 


2.37 


6.75 


0.37 


0.28 


P P7 [27] 


0.81 


0.53 


0.85 


2.19 


1.14 


1.99 


M pp + pj [28] 


0.29 


0.22 


0.35 


0.86 


2.86 


1.19 


MPR [29] 


0.25 


0.14 


0.24 


0.63 


4.0 


4.41 


SeSi [30] 


0.47 


0.26 


0.44 


1.18 


2.12 




5 Limits on 


the 


high 


energ 


;y 2), upper limits on the diffuse flux 



neutrino fluxes of extraterrestrial neutrinos are cal- 

culated. For a 90% confidence level 
an upper limit on the number of sig- 

Since no events have been observed nal events of ngQ% =2 . 5 is obtained 

which pass the final cuts (see Table 
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according to Conrad et al. [22] with 
the unified Feldman- Cousins order- 
ing [23]. We assume an uncertainty in 
signal detection of 24% and a back- 
ground of zero events (which leads to 
a conservative estimation of ngo% ac- 
cording to the Feldman-Cousins ap- 
proach) . 

The expected number of signal 
events N Vi for any assumed flux 
& Ui (E) of neutrinos of flavor i, is 
given by expression (5). A model of 
astrophysical neutrino sources, for 
which the total number of expected 
events, A r mo( ^ e |, is large than ngo%, 
is ruled out at 90% CL. iV mo( j e i is 
given as iV model = Y,N Vi . Table 4 
represents event rates and model re- 
jection factors (MRF) ^90%/-^model 
for models of astrophysical neutrino 
sources obtained from our search. 
Recently, similar results have been 
presented by the AMANDA collabo- 
ration [6,24], model rejection factors 
are shown in Table 4. 

The models by Stecker and Sala- 
mon [25] labeled "SS Q", as well as 
the models by Szabo and Protheroe 
[26] "SP u" and "SP 1" represent 
models for neutrino production in 
the central region of Active Galactic 
Nuclei. As can be seen from Table 
4, these models are ruled out with 

^90%/Nmodel ~ °- 06 " °- 4 - Further 
shown are models for neutrino pro- 
duction in AGN jets: calculations 
by Protheroe [27] and by Mannheim 
[28], which include neutrino produc- 
tion through pp and pj collisions 
(models "P prf and "M pp + prf , 
respectively), as well as an evalua- 
tion of the maximum flux due to a 
superposition of possible extragalac- 



tic sources by Mannheim, Protheroe 
and Rachen [29] (model "MPR" ) 
and a prediction for the diffuse flux 
from blazars by Semikoz and Sigl [30] 
"SeSi" . The latter models for blazars 
are currently not excluded. 

For an E~ 2 behaviour of the neutrino 
spectrum and a flavor ratio v e : : 
v T = 1 : 1 : 1, the 90% CL. up- 
per limit on the neutrino flux of all 
flavors obtained with the Baikal neu- 
trino telescope NT200 (1038 days) is: 

E 2 <$> < 8.1 x lO^cnrW^GeV. 

(8) 

For the resonant process 

v e + — > W~ — > anything, (9) 

with the resonant neutrino energy 
E = M 2 /2m e = 6.3 x 10 6 GeV and 
a cross section 5.02 x 10~ 31 cm 2 , the 
event number is given by: 

E+ 

N u - e =T J dVt J dEV eff J dE u x 

E~ 

x ^(E^N APH2 o^ (10) 

where E ± = (M w ± 2T w ) 2 /2m e and 
M w =80.22 GeV, Y w =2.08 GeV. 

Assuming an upper limit on the num- 
ber of signal events n 90 % =2.5, the 
model-independent limit on v e at the 
W - resonance energy is: 

^ < S.SxlO^cnrW^GeV" 1 . 

(11) 
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Fig. 12. Left panel: all-flavor neutrino flux predictions in different models of neu- 
trino sources compared to experimental upper limits to E~ 2 fluxes obtained by this 
analysis and other experiments (see text). Right panel: Baikal experimental limits 
compared to two model predictions. Dotted curves: predictions from model SS [25] 
and SeSi [30]. Full curves: upper limits to spectra of the same shape. Model SS is 
excluded (MRF=0.25), model SeSi is not (MRF=2.12). 



Figure 12 (left panel) shows our up- 
per limit on the all flavor E~ 2 diffuse 
flux (8) as well as the model inde- 
pendent limit on the resonant v e 
flux (diamond) (11). Also shown are 
the limits obtained by AMANDA 
and MACRO [6,24,31], theoreti- 
cal bounds obtained by Berezinsky 
(model independent (B) and for an 
E~ 2 shape of the neutrino spec- 
trum (B(E~ 2 )) [32], by Waxman and 
Bahcall (WB) [33], by Mannheim 
et al.(MPR) [29], predictions for 
neutrino fluxes from topological de- 
fects (TD) [30], prediction on dif- 
fuse flux from AGNs according to 
Nellen et al. (NMB) [34], as well as 
the atmospheric conventional neu- 
trino fluxes [35] from horizontal and 
vertical directions ( (z/) upper and 
lower curves, respectively) and at- 
mospheric prompt neutrino fluxes 
{y pr ) obtained by Volkova et al. [36]. 



The right panel of fig. 12 shows our 
upper limits (solid curves) on diffuse 
fluxes from AGNs shaped according 
to the model of Stecker and Salamon 
(SS) [25] and of Semikoz and Sigl 
(SeSi) [30], according to Table 4. 

The diffuse neutrino flux is assumed 
to be composed of contributions from 
cosmological sources with different 
luminosity and energy spectra. To 
benchmark the energy dependence 
of our limit without referring to a 
special model, we show in Fig. 13 the 
upper limits derived for E~ x spectra 
with variable, sharp cutoff. The parts 
of the spectra drawn with thick lines 
cover the energy ranges containing 
90% of expected events. Also shown 
is a curve connecting the points 
which correspond to the median en- 
ergies of recorded events. This curve 
represents the upper limit on a dif- 
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fuse flux formed by superposition of 
the benchmark spectra. Also shown 
(horizontal line) is the upper limit 
on an E~ 2 spectrum, as given in (8) 
above. 



-4 




3456789 10 



Ig(EyGeV) 

Fig. 13. Upper limits on E^ 1 bench- 
mark spectra with variable cutoff. Thick 
lines cover the energy ranges contain- 
ing 90% of expected events. Points at 
the median energies are connected by 
the thick curve. Also shown is the up- 
per limit on E~ 2 spectrum. 



6 Conclusion 

The neutrino telescope NT200 in 
Lake Baikal is taking data since 
April 1998. Due to high water trans- 
parency and low light scattering, the 
detection volume of NT200 for high 
energy u e , and v T events is several 
megatons and exceeds the geometri- 
cal volume by a factor of about 50 
for highest energies. This results in 
a high sensitivity to diffuse neutrino 
fluxes from extraterrestrial sources 
- more than an order of magnitude 
better than that of underground 



searches and similar to the published 
limits of AMANDA, the other oper- 
ating large neutrino telescope. The 
upper limit obtained for a diffuse 
{y e + + u T ) flux with E~ 2 shape is 
E 2 <$> = 8.1 x lO^cm-V^r^GeV. 
With 3.3 x lO-^cm-V^r^GeV" 1 , 
the limit on a i7 e flux at the resonant 
energy 6.3xl0 6 GeV is presently the 
most stringent. 

To extend the search for diffuse ex- 
traterrestrial neutrinos with higher 
sensitivity, NT200 was significantly 
upgraded. [37]. In March/April 2005 
we fenced a large part of the search 
volume (see Fig. 12, right) with three 
sparsely instrumented strings. The 
three-year sensitivity of this enlarged 
detector NT200+ on the neutrino 
flux of all flavors, with about 5 Mton 
enclosed volume, is approximately 
2 x lO^cm-W^GeV for E >10 2 
TeV, i.e. three-four times better than 
that of NT200. NT200+ will search 
for neutrinos from AGNs, GRBs and 
other extraterrestrial sources, neu- 
trinos from cosmic ray interactions 
in the Galaxy as well as high energy 
atmospheric muons with E^ > 10 
TeV. 
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